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We study the purely relaxational critical dynamics with non-conserved order 
parameter (model A critical dynamics) for three-dimensional magnets with dis- 
order in a form of the random anisotropy axis. For the random axis anisotropic 
distribution, the static asymptotic critical behaviour coincides with that of ran- 
dom site Ising systems. Therefore the asymptotic critical dynamics is governed 
by the dynamical exponent of the random Ising model. However, the disorder 
influences considerably the dynamical behaviour in the non-asymptotic regime. 
We perform a field-theoretical renormalization group analysis within the minimal 
subtraction scheme in two-loop approximation to investigate asymptotic and ef- 
fective critical dynamics of random anisotropy systems. The results demonstrate 
the non-monotonic behaviour of the dynamical effective critical exponent z e ff- 
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1. Introduction 



The concept of scaling plays a central role in modern theory of critical phenomena 
[P. Introducing a set of appropriate scaling variables a large amount of experimental 
and numerical data can be described by a few scaling functions [ EJ- The most 
prominent effect in dynamical critical phenomena is critical slowing down which 
consists in an increase of the relaxation time approaching the critical point. This is 
induced by the divergence of the correlation length £ at the critical point and causes 
also relaxation time r to diverge with the dynamical critical exponent z: 
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Renormalization group (RG) theory, with its concepts of invariance of the system 
at the critical point against changes in length scale gives a basis to universality in 
connection with fixed points and static [E] and dynamic scaling [0j. 

This holds however only in the asymptotic region in the vicinity of the critical 
point. Further away scaling breaks down and the description of critical phenomena 
becomes more complicated and involves non-universal characteristics both in statics 
and in dynamics [ d E] • 

Another important point for observing the behaviour in certain universality 
classes is the homogeneity of the system under consideration. Therefore the effect of 
impurities on critical behaviour is of considerable interest. However it turned out the 
disordered systems may show also scaling within a certain universality class. This 
might be the universality class of a pure system or a new one [ El El El E] • Moreover 
the changes introduced by the disorder depend on the type of this disorder; namely 
is it introduced by dilution (random site [ [TU] or random bond [ [TT] systems), or 
as a random field [ ^] , random connectivity [ El EI] or as an anisotropy [ • The 
defects may be correlated [0E] or not. It may even happen that the second order 
transition of the pure system is destroyed [ [TBI IT7j . 

It further turned out that considering a specific system with defects the behaviour 
near the critical point seems to be non-universal. Knowing that the non-asymptotic 
behaviour is non-universal it became necessary to study the non-asymptotic be- 
haviour of such systems in more detail. Indeed in many systems (e.g with site dis- 
order) effective critical behaviour can explain the experimental situation [ El d] • 

RG investigations of dynamic critical behaviour is in many cases technically much 
more involved in comparison to statics. As a consequence results for dynamics are 
known and in much lower approximations (in most cases only up to two loop order). 
On the other hand for the dynamics of magnetic systems with mode coupling terms 
the dynamical critical exponent is known exactly or contains only a static exponent. 

In this paper, we will present an analysis of the dynamical critical behaviour 
of random anisotropy magnets which constitute a large class of disordered systems 
[ I14j . In order to give some examples, the majority of the amorphous rare-earth con- 
taining alloys is recognized as random anisotropy magnets [ 13], certain crystalline 
compounds with an rear-earth component belong to this class too [H2j. Random 
anisotropy characterizes also molecular based magnets [ED], nanocrystalline materi- 
als [IS], as well as granular systems [ 122] . Moreover, the analysis involving random 
anisotropy found its application also in the interpretation of the phase transition in 
liquid crystals in porous media [E3]- 

The model currently used for the description of random anisotropy systems was 
introduced in the early 70-ies by Harris, Plischke, and Zuckermann [El] - It describes 
m-component spins located on the sites of a d- dimensional lattice, each spin sub- 
jected to a local anisotropy of random orientation. The Hamiltonian of the random 
anisotropy model (RAM) reads: 

U = - E JwSrS r , - D E^r^r) 2 . (2) 

R,R' R 

Here, Sr = (S' R , . . . , Sr), vectors R span sites of a <i-dimensional cubic lattice, 
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Dq > is an anisotropy constant, x is a random unit vector specifying direction 
of the local anisotropy axis. The interaction Jr.r/ is assumed to be ferromagnetic. 
Note, that for the Ising-like magnets, m — 1, the last term in (J2J) is just a constant, 
therefore the random anisotropy is present for m > 1 only. 

Below, we will consider quenched disorder, when the vectors £r, in (J2J) are ran- 
domly distributed with a distribution function p(x) and fixed in a certain configu- 
ration. It is well established by now, that the anisotropy axis distribution plays a 
crucial role for the origin of the low-temperature phase in the RAM. In particular, 
when the random vectors Xr, point with equal probability towards any direction, 
such a distribution may be called an isotropic one, the ferromagnetic ordering is im- 
possible [|2S] for spatial dimension d < 4. Whereas it may occur for an anisotropic 
distribution. In statics, this situation was corroborated by the RG studies of RAM 
[I2HJ 1211 I2H I2H EH] restricting x to point along one of the 2m directions of the axes 
ki of a hypercubic lattice (cubic distribution): 



with Kronecker deltas S(y). In this case, the second order phase transition into 
the magnetically ordered low-temperature phase occurs and in asymptotics it is 
characterized by the critical exponents of the random-site Ising model, the fact 
suggested already in Ref. [ 177] and confirmed later in Refs. [ [2E1 12H1 ED] • 

Taken that the studies of static criticality of random anisotropy magnets are far 
from being as intensive as those of the diluted magnets [ , even less is known 
about their dynamic critical behaviour. The dynamical models for systems with 
isotropic distribution of local anisotropy axis were considered in Refs. [ ETJ E21 EH] > 
the first order RG calculations were performed in Ref. [EI]- Whereas the problem of 
dynamic critical behaviour of a RAM with an anisotropic random axis distribution 
as to our knowledge up to now remained untouched. 

In this paper, we consider a purely relaxational dynamics of a three-dimensional 
(d = 3) RAM with non-conserved order parameter (model A in classification of 
Ref. [ E]) and a cubic random axis distribution ©. As far as the static critical 
behaviour of such magnets (note, with m > 1) belongs to the universality class 
of a random-site Ising model [ 123 123 EH EH] , for which the heat capacity does 
not diverge [CHI- Therefore, the critical dynamics of such a model is governed in 
asymptotics by the model A random-site Ising magnet dynamical critical exponent 
for any order parameter components number m. However in the non-asymptotic 
region, the model possesses a rich effective critical behaviour, as will be shown by 
our subsequent analysis. Taken that it is this effective behaviour which is observed 
both in experiments and in the MC simulations, it is important to have a RG 
prediction for typical scenarios of the approach to criticality in random anisotropy 
magnets. 

The rest of the paper is organized as following: in the next section |21 we present 
the Langevin equations governing model A dynamics and describe the renormaliza- 
tion procedure, in section El we give results of our calculations obtained in two- loop 




8=1 



(3) 
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approximation and display possible scenarios of effective critical behaviour. Conclu- 
sions and outlook are given in section |3J 



2. Model equations and renormalization 



We consider the dynamics for model (J2J) with random axis distribution (jSJ) to be 
relaxational with non-conserved m-component order parameter 0o = ipo(R). In this 
case the relaxation of the order parameter is described by the Langevin equation: 



dipt,. 



i,0 



1 . . . m, 



(4) 



dt d(p ifi 

with the Onsager coefficient T, stochastic forces 9 Vi obeying the Einstein relations: 
<6 Vi (R,t)6 ipj (R',t')> = 2V5(R-R')6(t-t')5 ij , (5) 
and the disorder- dependent equilibrium effective Hamiltonian TC: 



n 



d d R 



\V<p \ +r\(p \ J+^[l^o| -D{xtp 



(6) 



In (JEJ), the field <y?o is an m-component vector, D is proportional to the anisotropy 
constant of the spin Hamiltonian (J2J) with Dq] r and vq are defined by Dq and the 
fourth order coupling of the m- vector magnet (see Ref. [EH] for details). 

We treat the critical dynamics of the disordered model within the field theoretical 
RG method based on the Bausch-Janssen- Wagner formulation [OH]; where the ap- 
propriate Lagrangians are studied. For the model equations (jU)-© the Lagrangian 
reads: 



C = J d d RdtJ2&,o 



dt 5ip ifi 



(7) 



with a new auxiliary response m-component field <p Q . Here and below sums over field 
components span values from 1 to m. 

Studying critical properties of disordered systems one should average over the 
random degrees of freedom. In order to treat quenched disorder, in statics often the 
replica trick is used [G2]. However, it was established in [ 32j that in dynamics it is 
not necessary to make use of the replica trick: it is sufficient to the average over the 
random variables x with their distribution (J3J). Then the Lagrangian for the model 
reads: 



C 



d d RdtJ2&,o 

i 



+r(r -V 2 ) J (pi,o-Tcpi t Q+—(p i fl'^2 ( Pj,o ( Pj,o^ 



3! 



■<Pi 



+ 



(8) 



Here, tq is proportional to the temperature distance to the mean field critical point 
and the bare couplings are u > 0, v > 0, w < 0. Moreover, u and w are 
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connected to the moments of the distribution (jHJ) in such a way that wq/uq = —m. 
The term with x/q does not appear after averaging over disorder, however it should 
be added since it will be generated within the renormalization procedure. The set 
of static couplings {uq, t>o, Wo, yo} will be denoted below by {u 0ti }, i — 1, .., 4. 

G(k, uj)5(k + k')8{u + u/)6i d 
l,^ C(k, uj)8(k + k')5{u + ufjSij 



fc",l»" 




TA5(k + k' + k" + k"')5{uj + uj' + uj" + uj'") 
T 2 B5{k + k')5{k" + k'")8(u + uj')5(u" + uj'") 



Figure 1. Elements for construction of Feynman diagrams. The response function 
and the correlation function read: G(k,uj) = l/(—iuj + T(r + k 2 )), C{k,uj) = 
2r/| - iuj + T(r + A; 2 )| 2 . A stands for v /3! + 5^5 j %m + c>j,m<5j,i)/ 3 or for 

yo/3! Sij5j t i5i tm while 5 equals to u /3! 5i,j^ jm or to w /3\ Sijdj^i^- 

With the Lagrangian (jHJ) depending on the bare quantities (denoted by the sub- 
and superscripts V) we analyze within field theory [0] the dynamical vertex func- 
tions. As far as the static RG functions have been obtained before [E3 we need 
only to calculate the two-point dynamical vertex function r~^(r , {w^o}, f 1 , k, uj) = 
f &ip(ro, {wj,o}, f , k } The calculations are performed using Feynman diagrams, 

elements for them are given in the Fig. ^whereas the one- and two-loop contributions 
to T l £ are depicted in Fig. El 




Figure 2. Diagrams of the function f ^(ro, {uifl}, f , k, uj) up to two-loop order. 
First two terms represent one-loop contribution, while the rest of diagrams is of 
two- loop order. 



We perform renormalization of within the minimal subtraction scheme [Ej. 
In this scheme, in order to define renormalized static (p, r, {ui}) and dynamic (<£>, 
T) fields and couplings, the renormalization factors Z v , Z r , Z Ui and Z$, Z? are 
introduced by: 

<P = Z~ 1/2 f Q , ft = Z^ 1/2 ip , r = Z~ x r, Ui = Z' 1 Z 2 A d u , u T = Z r f. (9) 
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Here, \x is the external momenta scale, e — 4 — d, and Ad is a geometrical factor. 
The critical behaviour of the system is described by the following RG functions: 



MM) = 



CAM) 



-n- 



d\nZ r 



(r{M) = V 



d\nZ T 



(10) 



where the symbol means differentiation at fixed bare parameters. The /3-functions 
determine the RG flow of couplings under renormalization: 



^ = MM), 



1,..,4, 



and the flow parameter £ is related to the distance from the critical point. Subse- 
quently, an information about the critical behaviour of a system can be obtained 
from analysis of the fixed points (FPs) of the flow equations A FP {u*} is 

defined as simultaneous zero of all /3-functions: 



4. 



'121 



The stable and from initial conditions accessible FP corresponds to the critical point 
of the system. A FP is stable if all eigenvalues u>i of the stability matrix Bij = mis- 
calculated at this FP have positive real parts. 

The FP values of the RG ^-functions (jlOj) determine the asymptotic values of 
critical exponents. In particular, the dynamical asymptotic critical exponent z (pQ) 
is given at the stable and accessible FP by: 



z = 2 + ( r (u*,v*,w*,y*). 



(13) 



While the effective dynamical exponent z e s is calculated in the non-asymptotical 
region, where the renormalized couplings did not reach their FP values and it is 
defined by the solutions of the flow equations (fTTf : 

z eS = 2 + ( r (u{£),v{£),w(£),y{£)). (14) 

We neglect contributions to z c g coming from the amplitude function because they 
are considered to be small. 



3. Results 

The static RG functions within the minimal subtraction scheme are known in 
two-loop [E] approximation. Within the massive renormalization they have been 
calculated already in five-loop [E0| approximation. Calculating the dynamical func- 
tion £p within two- loop order we use the static RG functions of the same order [I15j. 
Since series for these static functions are known to be asymptotic at best we use 
Pade-Borel resummation scheme [ [32] in details described in Ref. [El- The FPs 
values and solutions of flow equations [ H2] obtained on the basis of these functions 
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we use in present study adding to them the new two-loop expression for the function 



Solving the FP equations f!12l) for the static /^-functions at fixed space dimension 
d = 3 [ 37] results in 16 FPs [ EI3 EH]. The region of physical importance u > 0, 
v > 0, w < includes 10 FPs. Below we list the most interesting FPs from them 
together with the asymptotic value for the z exponent (for the numerical values of 
the FPs coordinates obtained in two-loop approximation with the help of Pade-Borel 
resummation see [E3> the value of the z exponent, however, is calculated by a direct 
substitution of FP coordinates into Eq. (fE^ : 

• Gaussian FP I: u* = v* = w* = y* = 0; z{\/m) = 2; 

• pure FP II: v* ^ 0, u* = w* = y* = 0; z(m = 2) = 2.053, z(m = 3) = 2.051; 

• polymer FP III: u* ^ 0, v* = w* = y* = 0; z{\/m) = 1.815; 

• Ising FP V: y* ^ 0, u* = v* = w* = 0; z{Vm) = 2.052; 

• cubic FP VIII: u*^0, y*^0, u*=w*=0; z(m=2)=2.157, ^(m=3)=2.042; 

• Ising FP X: u* ^ 0, y* = -w*, v* = 0; z(Vm) = 2.052; 

• random Ising FP XV: w* ^ 0, y* = — w*, u* = v* = 0; z{\/m) = 2.139. 

Here, we keep the FP numbering of Refs. [ HH1 ESI EH EH ED] • From the above list, 
only the random Ising (XV) and polymer (III) FPs are stable. However the polymer 
FP is not accessible from the physical initial conditions. This leads to the conclu- 
sion [ E3 EM EH ED] that the random Ising FP XV governs the critical behaviour. 
Therefore, the m- vector magnets with cubic random axis distribution belong to the 
universality class of the random-site Ising magnets. The non-asymptotic critical be- 
haviour of the RAM differs essentially from that of the random-site model as was 
demonstrated in statics in Ref . [ EI] • The same concerns the non-asymptotic dy- 
namical critical behaviour: the critical slowing down in RAM is governed by z c g 
exponent as explained below. 

The crossovers between different FPs lead to a rich pictures of possible RG flows 
[ I15j. Many flows are influenced by the Ising FPs V and X. Introducing into (fT^j) 
several typical RG flows which start from the physical region of initial couplings one 
obtains different regimes for approach of the effective dynamical exponent z e fr to 
asymptotics. The dependence on the flow parameter i of z e g for easy-plane (m = 2) 
and Heisenberg (m = 3) magnets are shown in the Figs. El and 0] correspondingly. 
Flow 3 was chosen to be affected by both Ising FPs V and X, therefore both curves 
3 of Figs. El and E] demonstrate that a large region for z e s might exist with dynamical 



Cr- The last is derived from the vertex function T 



,l ~L discussed in Section El and reads: 



Cr 





(m + 2)uv + 5u 2 + 5w 2 + Wuw + 3yw + 3vw + Suy) . (15) 
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exponent values of the pure one-component (Ising) model A. Curves 6 correspond 
to flows which comes near the pure FP II and curve 7 to flows which come near the 
cubic FP VIII. 
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Figure 3. Effective critical exponent function of the logarithm of the 

flow parameter for order parameter dimension m = 2. Dashed line indicate the 
value of z at the FPs V, X. See text for details. 



Although the asymptotic exponents of the random anisotropy magnets consid- 
ered here are the same as those of the random-site (diluted) Ising magnets, the 
approach to the asymptotical region essentially differs from the diluted magnets. It 
is defined by smallest static stability exponent uj = —0.0036 [ which is equal 
in absolute value to the ratio of heat capacity critical exponent a r and correla- 
tions length critical exponent u r of random-site Ising model [EDI- As a consequence 
the Wegner correction to scaling is A = uu r = —a r . The high-loop estimate gives 
A m 0.049 ± 0.009 [|3U]. Such a small value of A means that the approach to the 
asymptotic values is very slow. Therefore practically only the non-asymptotic critical 
behaviour governed by effective critical exponents will be observed experimentally 
or in the numerical simulations. 

As it is seen from Figs. El IH another particular feature of z e s seems to be that 
it reaches its asymptotic value z always from the region z e s > z. Therefore, in an 
experimental situation, a decrease of z e s may serve as an evidence of approach to the 
asymptotics. Note that such a scenario is an intrinsic feature of the critical slowing 
down in random-anisotropy magnets. When disorder is implemented by dilution 
of the non- magnetic component, an approach of z e R to its asymptotic value is not 
necessarily only from above [I38j. 



4. Conclusions 



In this paper, we have analyzed the critical slowing down in magnets influenced 
by random anisotropy. These magnets have a second order phase transition to ferro- 
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In / 



Figure 4. Effective critical exponent ZqQ ELS cl function of the logarithm of the 
flow parameter for order parameter dimension m = 3. Dashed line indicate the 
value of z at the FPs V, X. See text for details. 

magnetic order for an anisotropic (cubic) random axis distribution. Therefore, our 
goal was to study relaxational dynamics of the non-conserved order parameter in the 
vicinity of the phase transition point. For this purpose we completed previous static 
RG calculations [[12112211201 by calculating the two- loop dynamical RG function £ r 
given in Eq. (|15j) . Combining this result with the former data for the static critical 
behaviour we obtained numerical values for the effective critical exponent z e s which 
governs the critical slowing down of the relaxational time when T c is approached. In 
Figs. 13141 we give results for two most physically interesting cases m = 2 and m — 3, 
which correspond to the easy-plane and Heisenberg random anisotropy magnets. 

Although the asymptotic dynamical critical behaviour of random anisotropy sys- 
tems with cubic distribution is the same as for the random-site Ising systems, the 
crossover between different fixed points considerably influences the non-asymptotic 
critical properties. Different scenarios of dynamical critical behaviour are observed. 
Since the approach to asymptotics is very slow, it might be observed in real and nu- 
merical experiments. The effective exponents measured may take values essentially 
differing form the asymptotic one (in our calculation z = 2.139). For example in a 
large region z e g can be equal to the exponent of the pure Ising model (z = 2.052). 
One more particular feature of critical slowing down in random anisotropy magnets 
which is predicted by our analysis is that, contrary to the diluted magnets, z e fr seems 
to reach its asymptotic value z always from the region z c s > z. 

Another important contribution to the effective dynamical exponent could come 
from the coupling of the order parameter to a conserved density (changing from 
model A dynamics to model C [Ej). Since the stable fixed point has a non-diverging 
specific heat the asymptotic discussed here would not be changed [ EH] ■ A more 
detailed account on that is in preparation. 

It is our pleasure and honour to contribute this paper to the conference on the 
occasion of Prof. I.R. Yukhnovskii 80th birthday. His early work on the phase tran- 
sitions theory and on an account of the non-asymptotic criticality preceded in 
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many respects many later contributions to this field. One of us (Yu.H.) is deeply 
indebted to the jubilee for introducing into the fascinating field of phase transi- 
tions and critical phenomena. R. F. acknowledges the fruitful cooperation with the 
Institute for Condensed Matter Physics. 

This work was supported by Austrian Fonds zur Forderung der wissenschaftli- 
chen Forschung under Project No. PI 65 74. 

References 

1. C. Domb. The critical point. Taylor & Francis, London, 1996; P.M. Chaikin, T.C. 
Lubensky. Principles of condensed matter physics. Cambridge University Press, Cam- 
bridge, 1995. 

2. B. Widom. J. Chem. Phys., 1965, 43, 3898; L.P. Kadanoff. Physics, 1966, 2, 263; A.Z. 
Patashinskii, V.L. Pokrovskii, Zh. Eksp. Teor. Fiz., 1966, 50, 439 [Sov. Phys. JETP, 
1966, 23, 292]. 

3. J. Zinn-Justin. Quantum Field Theory and Critical Phenomena. Oxford University 
Press, Oxford, 1996; H. Kleinert, V. Schulte-Frohlinde. Critical Properties of c/> 4 - 
Theories. World Scientific, Singapore, 2001. 

4. B. I. Halperin, P. C. Hohenberg, Rev. Mod. Phys., 1977, 49, 436. 

5. I.R. Yukhnovskii, Theor. Math. Phys., 1978, 36, 798 [Teor. Mat. Fiz., 1978, 36, 373]; 
I.R. Yukhnovskii. Phase transitions of the second order. Collective variables method. 
World Scientific, Singapore, 1987. 

6. A.B. Harris, J. Phys. C, 1974, 7, 1671. 

7. D. Boyanovsky, J.L. Cardy, Phys. Rev. B, 1982, 26, 154; Phys. Rev. B, 1983, 27, 6971. 

8. A. Weinrib, B.I. Halperin, Phys. Rev. B, 1983 27, 413. 

9. J.M. Luck, Europhys. Lett., 1993, 24, 359. 

10. For the recent reviews on criticality of diluted magnets see e.g.: A. Pelissetto, E. Vicari, 
Phys. Rep. 2002, 368, 549; R. Folk, Yu. Holovatch, T. Yavors'kii, Physics - Uspekhi, 
2003, 46 169 [Uspekhi Fizicheskikh Nauk, 2003, 173, 175]. 

11. P. E. Berche, C.Chatelain, B. Berche, W. Janke, Eur. Phys. J. B, 2004, 39, 463; B. 
Berche, P. E. Berche, C.Chatelain, W. Janke, Condens. Matter Phys., 2005, 8, 47. 

12. D. P. Belanger, A. P. Young, Jour. Magn. Magn. Mater., 1991, 100, 272. 

13. W. Janke, M. Weigel, Phys. Rev. B, 2004, 69, 144208. 

14. A review of early work on random anisotropy magnets may be found in: R. W. 
Cochrane, R. Harris, M. J. Zuckermann, Phys. Reports, 1978, 48, 1. Recent experi- 
mental, numerical, and theoretical studies are reviewed in Ref. [I15j. 

15. M. Dudka, R. Folk, Yu. Holovatch, Journ. Mag. Mag. Mat., 2005, 294, 305. 

16. Y. Imry, S.-k. Ma, Phys. Rev. Lett., 1975, 35, 1399. 

17. R.A. Pelcovits, E. Pytte, J. Rudnick, Phys. Rev. Lett., 1978, 40, 476; S.-k. Ma, J. 
Rudnick, Phys. Rev. Lett., 1978, 40, 589. 

18. M. Dudka, R. Folk, Yu. Holovatch, and D. Ivaneiko, J. Magn. Magn. Mater., 2003, 
53, 243, A. Perumal, V. Srinivas, V. V. Rao, R. A. Dunlap, Phys. Rev. Lett., 2003, 
91, 137202; R. Folk, Yu. Holovatch, T. Yavors'kii, Phys. Rev. B, 2000, 61 15114, P. 
Calabrese, P. Parruccini, A. Pelissetto, E. Vicari, Phys. Rev. E, 2004, 69, 036120. 

19. A. del Moral, C. de la Fuente, J. I. Arnaudas, Phys. Rev. B, 1996, 54, 12245; J. Phys., 
1996, 8, 6945. 



10 



Critical slowing down in random anisotropy magnets 



20. M. A. Girtu, C. M. Wynn, J. Zhang, J. S. Miller, A. J. Epstein, Phys. Rev. B, 2000, 
61, 492. 

21. P. Marin, M. Vaazquez, J. Areas, A. Hernando, J. Magn. Magn. Mater. 1999, 203, 6; 
M. Zheng, R. Skomski, Y. Liu, D. J. Sellmyer, J. Phys., 2000, 12, L497; M. L. Bue, 
V Basso, C. Beatrice, P. Tibero, G. Bertotti, J. Magn. Magn. Mater., 2002, 242-245, 
1089. 

22. W. C. Nunes, M. A. Novak, M. Knobel, and A. Hernando, J. Magn. Magn. Mater., 
2001, 226-230, 1856. 

23. J. Chakrabaty, Phys. Rev. Lett., 1998, 81, 385. 

24. R. Harris, M. Plischke, M. J. Zuckermann, Phys. Rev. Lett., 1973, 31, 160. 

25. As shown in the Ref. [EI; an absence of the ferromagnetic ordering for an isotropic 
random axis distribution at d < 4 follows from the Imry-Ma arguments first formulated 
in the context of the random field Ising model in Ref. [ I16j . Later this fact has been 
proven by several other methods (see [ 29 for a more detailed account). 

26. A. Aharony, Phys. Rev. B, 1975, 12 1038. 

27. D. Mukamel, G. Grinstein, Phys. Rev. B, 1982, 25 381. 

28. A. L. Korzhenevskii, A. A. Luzhkov, Sov. Phys. JETP, 1988, 67, 1229 [Zh. Eksp. Teor. 
Fiz., 1988, 94, 250]. 

29. M. Dudka, R. Folk, Yu. Holovatch, Condens. Matter Phys., 2001, 4, 77; ibid., 2001, 
4, 459. 

30. P. Calabrese, A. Pelissetto, E. Vicari, Phys. Rev. E, 2004, 70, 036104. 

31. S.-k. Ma, J. Rudnick, Phys. Rev. Lett., 1978, 40, 587. 

32. C. De Dominicis, Phys. Rev. B, 1978, 18, 4913. 

33. A. Khurana, Phys. Rev. B, 1982, 25, 452. 

34. U. Krey, Z. Phys. B, 1977, 26, 355. 

35. R. Bausch, H. K. Janssen, and H. Wagner, Z. Phys. B, 1976 24, 113. 

36. G. A. Baker, Jr., B. G. Nickel, D. I. Meiron, Phys. Rev. B, 1978, 17, 1365. 

37. R. Schloms, V. Dohm, Europhys. Lett., 1987, 3, 413; Nucl. Phys. B, 1989, 328, 639. 

38. V. Blavats'ka, M. Dudka, R. Folk, Yu. Holovatch. to appear in Phys. Rev. B, 2005; 
preprint cond-mat/0506644 

39. M. Dud ka, R. Folk, Yu. Holovatch, G. Moser, submitted to Phys. Rev. E, 2005; preprint 
cond-mat/0506325| 



11 



